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Abstract. We discuss the subtle relationship between massive gravity and bimetric 
gravity, focusing particularly on the manner in which massive gravity may be viewed 
as a suitable limit of bimetric gravity. The limiting procedure is more delicate than 
currently appreciated. Specifically, this limiting procedure should not unnecessarily 
constrain the background metric, which must be externally specified by the theory 
of massive gravity itself. The fact that in bimetric theories one always has two sets 
of metric equations of motion continues to have an effect even in the massive gravity 
limit, leading to additional constraints besides the one set of equations of motion 
naively expected. Thus, since solutions of bimetric gravity in the limit of vanishing 
kinetic term are also solutions of massive gravity, but the contrary statement is not 
necessarily true, there is not complete continuity in the parameter space of the theory. 
In particular, we study the massive cosmological solutions which are continuous in the 
parameter space, showing that many interesting cosmologies belong to this class. 
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1. Introduction 

Massive gravity has recently undergone a significant surge of renewed interest. Since 
de Rham and Gabadadze [I], (see also de Rham, Gabadadze, and Tolley [2]), 
demonstrated that it is possible to develop an extension of the Fierz-Pauli mass term for 
linearized gravity [3], one that avoids the appearance of the Boulware-Deser ghost [I], 
at least up to fourth order in non-linearities, activity on this topic has become intense 
— with over 50 articles appearing in the last two years. 

Subsequently, both Hassan and Rosen [5], and Koyama, Niz, and Tasinato j6], 
have independently re-expressed and re-derived the theory of de Rham and Gabadadze 
- simplifying the treatment and shedding additional light on its characteristics. In 
fact, Hassan and Rosen took a significant step further by extending the theory to a 
general background metric [5]. Rapidly thereafter, several papers appeared studying 
the foundations of this theory [7J [9J [101 HH [12], and proving the absence of ghosts 
in the nonlinear theory [13j El EE]. In particular, Hassan and Rosen showed that 
the introduction of a kinetic term for the background metric in the ghost-free massive 
gravity leads to a bimetric gravity theory which is also ghost-free [17J. It should be noted 
that in this case the background metric is not only an externally-specified kinematic 
quantity, but also has its own dynamics, acquiring the same physical status as the 
"foreground" metric [18]. More recently, ghost-free multi-metric theories have also been 
considered [T9| 120]. 

Any cautious and physically compelling approach to massive gravity should not only 
respect the beauty of standard general relativity, but also enhance it in some manner, by 
embedding standard general relativity in some wider parameter space. Physically, one 
might hopefully expect that the observational predictions of the extended theory should 
be continuous in these extra parameters, and that general relativity would be recovered 
by taking the limit for a zero graviton mass in massive gravity theories. However, this is 
not necessarily the case, since (as is well known) the predictions of massive gravity often 
qualitatively differ from those of general relativity even when the graviton mass vanishes. 
This effect is known as the van Dam-Veltman-Zakharov (vDVZ) discontinuity |21|. [22] . 
The Vainshtein mechanism provides us (under some appropriate conditions) with a 
way to avoid (or rather ameliorate) such a discontinuity [23] (see also [2H ESI [251 [27] 
and [6]). Furthermore, as has been pointed out in reference [5], the vDVZ discontinuity 
can sometimes be arranged to be absent when one considers the background metric to 
be non-flat in massive gravity, which is certainly the more general situation. That is, 
it seems that the predictions of massive gravity might be arranged to be continuous in 
parameter space when suitably taking into account a curved background [28, 29J — this 
is one of several reasons for being interested in non-trivial backgrounds. 

But why, in the first place, should one modify general relativity in such a manner 
and even entertain the possibility of massive gravity? References [28] and [29] provide 
historical reviews of the motivations. On the other hand, over the last decade the 
theoretical revolution based on the inferred accelerated expansion of the universe has 
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become a good reason for considering a wide variety of possible modifications of general 
relativity [301 EH [33J |32j EI]- Consequently, when de Rham, Gabadadze, and Tolley 
found a family of ghost-free (flat-background) massive gravity theories [U |2], their 
cosmological consequences were quickly analyzed [351 ESI EH EH], even though the 
theory had not yet been formulated to be compatible with a non-flat background 
metric. Although these studies can be considered as a first attempt to understand the 
cosmological consequences of ghost-free massive gravity, their results are severely limited 
by the (with hindsight unnecessary, and perhaps even physically inappropriate) choice 
of a flat background metric. Indeed, assuming a flat background is not the most general 
situation, and it is even contra-indicated when considering black hole geometries [13] or 
cosmological scenarios [39J. Later on, independently and almost simultaneously, three 
groups have considered cosmological solutions in bimetric gravity [40], UH [42] , while also 
studying massive-gravity cosmologies more or less in passing. However, as we will show 
in this paper, considering solutions of massive gravity as a limit of bigravity theories 
can be much more subtle than expected. 

At this point one might still reasonably wonder whether bimetric gravity is itself 
well-behaved in parameter space. Although a background kinematic metric has provided 
us with a continuous limit of massive gravity with respect to general relativity, it could 
still be that the introduction of a kinetic term for the background metric, (with its 
associated dynamics), disrupts coherence with respect to massive gravity. In other 
words, it may well be that the physical predictions of bimetric gravity differ from those 
of massive gravity in the limit where this theory should be recovered. To settle this issue, 
in this paper we study the massive gravity limit of bimetric gravity in full generality. 
As we will show, the solutions of bimetric gravity in the limit where the kinetic term 
for the background metric vanishes will also be solutions of massive gravity compatible 
with a non-flat background metric, but not necessarily vice versa. 

This paper is organized as follows: In section [2] we (briefly) summarize some 
previous results of massive gravity and bimetric gravity. In section El we explore 
various limits by which bimetric gravity might be used to reproduce massive gravity. 
In particular we investigate how the vanishing kinetic term limit should be considered, 
see section 13.11 concluding that the limit procedure is not implying the need for a flat 
background metric in massive gravity. Consideration of this limit allows us to obtain the 
cosmological solutions which are continuous in the parameter space of massive gravity, 
see section 13.21 We study these solutions which can be classified as general solutions, 
continuous cosmological solutions of any ghost-free theory, and special case solutions. 
Although the second group are solutions only for a particular kind of massive gravity 
models, we show that they present some features of particular interest. In section [3731 
we briefly comment the consequences of retaining some effects of the background matter 
when taking this limit. Section H] contains some discussion and conclusions, while some 



purely technical formulae and computations are relegated to Appendix A[ |Appendix B 



and Appendix C 
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2. Massive gravity and bimetric gravity 

As is already rather well-known [3] (see also |39] ) . in order to consider massive gravity 
one is forced to introduce a new rank-two tensor whose kinematics is at this stage 
externally specified, and not governed by the theory itself. This new rank-two tensor f^ v 
can best be interpreted as a "background" metric, not necessarily flat, with linearized 
fluctuations h^ u of the physical "foreground" metric g^ u = f^ u + h^ u satisfying a massive 
Klein-Gordon equation [4j (see also [39] and [33] ) - That is, the (full non-linear) mass 
term appears in the action as some interaction term which depends algebraically on the 
tensors f^ u and g^ u , through the quantity 

(g-'fYv = g^fpu, (i) 

and which in the linearized limit reproduces a suitable quadratic mass term. On the 
one hand, there is a non-denumerable infinity of such interaction Lagrangians, on the 
other hand, almost all of them lead to the physically unacceptable Boulware-Deser 
ghost [4]. For example, a particular specification of the interaction Lagrangian given in 
reference [39] led to the cosmological models presented in references [33 SHI S3 EE], 
which assumes a flat background metric contraindicated for this situation, and 
reference [50] , which additionally contains specific technical criticisms to (33 SHJ H3 HE] . 
Different approaches include the Lorentz-violating massive gravity [4"9| |5T| 152] and the 
3-dimensional "new massive gravity" [531 E3] • The novelty of the last two years is that 
the interaction Lagrangian has now been very tightly constrained and fixed to depend 
only on a finite number of parameters |5| by requiring the theory to be ghost-free. 
A quite unexpected technical result (from the Boulware-Deser point of view [3], see 
also [39]) is that the dependence of this ghost-free interaction term on the background 
and foreground metrics is through the square-root quantity 

= (VivYu, that is 7^7% = g^U- (2) 

On the other hand, bimetric gravity was first introduced by Isham, Salam and 
Strathdee [18] (see also [55], EH EZ]) to account for some features of strong interactions, 
and it was later rejuvenated by Damour and Kogan in order to address new physics 
scenarios [58] (see also [591 EQl El]). It has also been recently proven to be ghost-free 
when considering the same interaction term as for massive gravity [T7] . 

Let us now focus our attention on 4-dimensional (and Lorentz invariant) massive 
gravity. We can express the action generally as 

Smg = ~T^G J /a; v /Z 5{%) + 2A-2 m 2 L mt ( 9 - 1 /)} + 5 w , (3) 

with S(xa) describing the usual matter action, with matter fields coupled only to the 
foreground metric g^, (and the measure y/—g ), to agree with the Einstein equivalence 
principle. The parameter m sets the scale for the graviton mass, and the interaction 
term L int ((7 -1 /) is a scalar chosen to be dimensionless. 

Regarding bimetric gravity, in addition to the kinetic term of the background 
metric, one must also consider the possibility of a background cosmological constant 
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A, and "background matter" Sr m \ coupling to f^ u . This now leads to 

Sbg = - J d 4 x^ {R(g) + 2 A - 2m 2 L int } + S M 

/^V /Z 7{W) + 2A}+^ (m) . (4) 

The effective Newton constant for the background spacetime/matter interaction is eG/n. 
The two parameters k and e can in principle be adjusted independently. 

The interaction term of the ghost-free theories can without loss of generality be 
written as [SJ [T7| (see also the discussion of Appendix B ) 

Lmt = e 2 (K) - c 3 e 3 (K) - c 4 e A (K), (5) 

with 

K^ = 5" u -Yu, (6) 



and the polynomials (see Appendix A for a more formal definition and properties) 
are 

e2(K) = ~ ([K] 2 - [K 2 ]) ; (7a) 
e 3 (K)= 1 -([K] 3 -3[K][K 2 ]+2[K 3 ]) ] (7b) 

UK) = 1 {[Kf - 6[K 2 ][K] 2 + 3[^ 2 ] 2 + 8[K][K*] - 6[iv 4 ]) ; (7c) 

where [K] = ti(K^ v ), and our definition of K, equation ([6]), agrees with that of 
references [B] and [40J, but differs from that of reference [5]. One can now consider 
the foreground tetrad e^, and the background inverse tetrad w^ A defined by [40] (see 
also references ESI El HH]) 

<T = ^ B e^ A e v B , fju, = Vab w„ A w v B , (8) 

where no direct analogy between these expressions and the definition of the Stiickelberg 
fields should necessarily be deduced (the relation between quantities defined in the 
tangent space and in the tetrad basis cannot be thought of as in any way recovering any 
gauge freedom). This formalism allows one to write the term involving the square root 
as 

-f v = eP A w\, (9) 

by requiring the consistency condition 

e^A w Bf > = e^s w Afl , (10) 

which leaves the equations of motion unchanged [40J. 

The equations of motion for massive gravity are obtained by varying the action ([3]) 
with respect to g^ v . These are 

G* v - A 5» v = m 2 T (efT) ^ + 8ttG T (m) ^, (11) 
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with T( m ^ u denoting the usual stress-energy tensor of the matter fields, while 

T^ u = r^-S^L int , (12) 

is the dimensionless graviton-mass-induced contribution to the stress-energy. It can be 
noted that t^ v can be written as 

t m -v* 9Lint - c»„ 9Lint (11) 

T p dr P ~ B de B ' {16} 

for any interaction Lagrangian which depends on the metrics through terms of the form 
[K n ], with n an integer, (or equivalently [7™]). In addition, the Einstein tensor of 
equation (ITT]) needs to satisfy the (contracted) Bianchi identity Taking into account 
the invariance of Sr m -\ under diffeomorphisms, which implies V^T*™^ = 0, that identity 
leads to a constraint on the graviton-mass-induced effective stress-energy: 

V M T (eff) ^ = 0. (14) 

On the other hand, the bimetric gravity theory now has two sets of equations of 
motion, obtained by varying with respect to the two metrics. Thus, in addition to 
equation ( TTTT) one has 

k (0"„ - A 5%) = m 2 T\ + e 8irG f (m) ^, (15) 

where the effective stress-energy tensor for the background metric is [ffij H2] 

r\ = t\. (i6) 

Furthermore, the Bianchi-inspired constraint which follows from equation f jT5|) is 
equivalent to that already obtained in equation ( IT4"]) . see [4T)| I4T~]. 

It should be emphasized that the equations of motion (ITT]) . ( TT5|) the definition of the 
effective energy-momentum tensor, equations (|T2|) . (|T3|) . (|T6|) and the constraint (EH]), 
are all completely independent of the particular form of the interaction term, this being 
an automatic result of the fact that the interaction term depends only on the quantity 
g~ x f through terms of the form [7™]. See for instance reference [39]. (Warning: There is 
a significant typo in reference [39] , amounting to accidentally dropping the scalar factor 
\[JJg in the effective stress-energy due to the graviton mass term — see reference [50J - 
fortunately this does not quantitatively affect weak-field physics, nor does it qualitatively 
affect strong-field physics — though it will certainly significantly change many of the 
details.) 

If we now consider the specific family of ghost-free theories given by ([5]), we can 
explicitly express t^ v in terms of matrices as 

r = ([ 7 ] - 3)7 - 7 2 + c 3 (e 2 (^)7 - ei(K) 7 • K + 7 • K 2 ) 

+c 4 (e 3 (A-)7 - e 2 {K) 1 ■ K + e^K)-/ ■ K 2 - 7 • K 3 ) . (17) 

An equivalent index-based formula can be found in reference [4"0] . 
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3. Continuity of massive gravity with respect to bimetric gravity 

An implication of the discussion in the previous section is that if one considers some 
particular solutions of a bimetric gravity theory, then those solutions will be more 
constrained than those corresponding to a massive gravity theory with the same 
interaction term. That is because when one considers f^ v to be non-dynamical (being 
externally specified by the definition of the theory), the kinetic term for this metric is 
no longer present and one should not consider the variation of the action with respect 
to this metric. 

Following this spirit, von Strauss et at have considered in reference [H] cosmological 
solutions for the ghost-free bimetric and massive gravity theories, given by the actions (111) 
and ([3]), respectively, and the interaction term ([5]). They have noted consistently that in 
massive gravity one loses the equations of motion given by equation (IT5l) . Nevertheless, 
in references [10| W2\ the authors also study cosmological solutions of massive gravity but 
considering it as a particular limit of bigravity theory. They have obtained different and 
non-equivalent results to that presented in jjl] because when considering massive gravity 
as a limit of bimetric gravity, one has not only the equations of motion of the foreground 
metric but also additional constraints. Thus, one may reasonably wonder whether the 
consideration of this limit implies that there is some kind of physical discontinuity in 
the parameter space of massive gravity. 

On the other hand, if one wants to consider massive gravity as a limit of bimetric 
gravity, then this limit should be carefully taken to avoid inconsistent results. In several 
references, presenting both pre-ghost-free and ghost-free analyses, the authors have 
concluded or implied that the background metric should be flat (Riemann-flat), an 
interpretation that can hide some problems of the theory, since a flat metric is an 
incompatible background for the cosmological scenarios that they have studied. The 
recent models for ghost-free massive gravity considered in references [lOl H2] amount in 
the current language (and now including the possibility of a background cosmological 
constant) to taking the limit k — > oo while holding m fixed and setting e = 0. The 
background equations of motion then degenerate to Q^ u — A/ Mi , = 0, so that the 
background spacetime is some Einstein spacetime — for example, Schwarzschild/ Kerr/ 
de Sitter/ anti-de Sitter, or even Milne spacetime, not necessarily Minkowski spacetime. 
(Implicitly ignoring any background cosmological constant, as in jlQl H2], leads to a 
Ricci-flat space for finite 7^,.) However, following the philosophy of massive gravity, the 
background metric should not be constrained by the limiting procedure, but externally 
specified by the theory itself. Thus, if one wants to recover massive gravity as a limit 
of bimetric gravity, then it seems inconsistent to consider a limit which (unnecessarily) 
fixes the background metric. 

3.1. The non- dynamical background limit 

If we consider the action of bimetric gravity given by equation (J4]), it is easy to see that 
one can recover the action of massive gravity, equation ([3]), simply by simultaneously 
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taking the limits k — > and e — > 0. Furthermore, the consideration of the limit «; — > 
will imply some constraints on the kind of interaction between both gravitational sectors 
in bimetric gravity, leading to one that must be also compatible with massive gravity to 
recover this theory, while not fixing the background metric as in references [IHl H2], at 
least in principle. For consistency it can be checked that the same constraints coming 
from the consideration of this limit can be obtained by taking directly the variation of 
the action of massive gravity ([3]) with respect to the background metric f^ u . Therefore, 
if one were interested in recovering the predictions of massive gravity in a certain limit 
of bigravity theory, then it would be natural to consider such a limit, which corresponds 
to the limit of vanishing kinetic term (that is the limit of a vanishing effective Planck 
mass for one of the metrics) . 

Considering k — > and e — > in the background equation of motion ( fl5l) . and 
using the definition of embodied in equation ( Fl6i) . the following perhaps unexpected 
constraint can be obtained: 

t\ = 0. (18) 

Moreover, taking into account equation ffl2l) . one can note that equations ffl8|) imply 
that the effective energy-momentum tensor appearing in the equations of motion for the 
dynamical metric, equation ( fill) , can be written as 

T^K = -8» v L int . (19) 

Thus, the constraint ( JT4l) . ultimately coming from the contracted Bianchi identity, 
implies 

<9 A Lint = 0. (20) 

Therefore, from equations (1T91) and ( 1201) we can conclude that the modification in the 
equations of motion (flT!) due to a putative non zero-mass for the graviton, considered 
as a limit of a theory with two dynamical metrics, is equivalent from the point of view 
of the physical foreground metric to simply introducing a cosmological constant, one 
which can at least in principle be either positive or negative. That is 

= -5^ A eff . (21) 

On the other hand, although in the next section we will focus our attention on theories 
of massive gravity with a spherically symmetric background metric, it should be pointed 
out that we have obtained no restriction about the curvature or dynamics of the 
background metric. In particular, there is no requirement for this metric to be Riemann 
flat (or Ricci-flat, or even an Einstein spacetime). It could be argued that one can choose 
to consider massive gravity with an Einstein background metric, and one can in fact 
do it. Nevertheless, that cannot be a requirement coming from seeing this theory as a 
particular limit of bimetric gravity, because in that case the adoption of that particular 
point of view would change the philosophy of the theory itself. 

Finally, it should be noted that up to now, we have not assumed any particular 
form for the interaction term (only that it depends on the metrics through [7"]), neither 
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have we assumed any symmetry for the metrics. Therefore, these results are completely 
general when studying the solutions for any massive gravity theory of that class as a limit 
of bimetric gravity. Thus, one can already conclude that there are fewer solutions that 
are solutions simultaneously of massive gravity and of bimetric gravity in the limit k — > 
and e — > 0, than of massive gravity alone. This is because the former solutions must be 
also solutions of massive gravity and in addition fulfill some extra constraints, namely 
equations (fT8l) . (fT9l and fl20l) . Moreover, those solutions, if any, would be equivalent to 
simply considering a foreground cosmological constant. 



3.2. Ghost-free cosmologies 

We now focus our attention on the ghost-free case, which corresponds to the specific 
interaction term given by equation ([5]), and study solutions of massive gravity which 
are continuous in the parameter space (that is solutions of both massive gravity and 
bimetric gravity in the non-dynamical background limit). As we have concluded that 
those solutions mimic the effect of a cosmological constant in the foreground space, it is 
of particular interest to consider cosmological solutions. Thus, we assume a spherically 
symmetric situation, where both metrics fulfill this symmetry. Both for massive gravity 
and for bigravity theories these metrics can be written in general as [40] 



/ M „ dx" dx u = {Adt + Cdrf - [Bdr- ^dt ) - U 2 dn 2 {2) , (23) 



2 



g„ u dx» dx v = S 2 dt 2 - N 2 dr 2 - R 2 dVl 2 2) , (22) 

and 

~N 

where all the metric coefficients are functions of t and r. 

As is well known [56] \57\ 160] (see also reference [40] for the ghost-free theory) , in 
bimetric gravity the requirement of T^ eS \ r = (or, equivalently, r° r = 0), which is 
implied by the consideration of a spherically symmetric scenario, leads in general to two 
classes of solutions: (i) those in which both metrics can be written in a diagonal way 
using the same coordinate patch, and (ii) those with metrics which are not commonly 
diagonal. Nevertheless, when one considers solutions of bimetric gravity with a vanishing 
kinetic term, such solutions must fulfill constraints f lTS]) instead of the equations of 
motion of the background metric (TT5T) . leading Q^ v unspecified and, therefore, being 
compatible with massive gravity. This requirement leads to the conclusion that there 



is no Lorentzian-signature solution for equations (1181) if C ^ (see Appendix C for 
details) without conflicting with [56j E3 ES HQ] • So we can simplify the background 
metric and take 

/ M „ dx^ dx v = A 2 dt 2 - B 2 dr 2 - U 2 dtt 2 {2) , (24) 

Considering metrics (122]) and (124]) . it can be seen that there are two general solutions, 
and additionally two special-case solutions, to equations (ITSl) and (|2"0"]) for any ghost-free 
massive gravity theory whose interaction term can be described by ([5]). 
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3.2.1. First general solution. The first general solution we discuss is particularly 
interesting, since it relates the two metrics in a very simple way - - a position- 
independent rescaling. The metrics satisfy 

/ M „ dx" dx v = D 2 g^ dx^ dx u , (25) 

where D is a constant with D ^ 1 such that 

C4 " {D-iy ■ (26) 

Explicitly 




D = l + — ±±W l + — ± -1. (27) 



For this solution the graviton mass produces a term in the modified Einstein equations 
(TTTj) . which can be described by the effective stress tensor 

T (eff)^ = _ [3 + C3 ( D _ !)] ( D _ !)2 ^ (28) 

which mimics the behavior of a positive cosmological constant if 03(1 — D) > 3, and 
a negative one otherwise. Therefore, these solutions can describe a universe which 
is accelerating as if this acceleration would be originated by a cosmological constant. 
Moreover, this foreground universe would have exactly the same symmetry as the 
background metric — that is having a homogeneous and isotropic foreground universe 
is possible only if the background metric is also FLRW with the same sign of spatial 
curvature. 

We might have argued on general grounds for the existence of a solution of this type 
for a generic interaction term and in the absence of any symmetry of the metrics once 
taking into account equation f|T8|) . When calculating t^ u , notice that it is proportional 
to dLi nt /dj u M . Thus, one might guess the existence of some general (not necessarily 
unique) solution with 7^ = D 8^ v , where the value of D is fixed by the theory. In fact, 
in view of the specific form of L int , one might reasonably infer a generic polynomial 
constraint on D. (That is because r cx dL- mt /d^ ~ Yli^ e il^l ~ ^kVk{.l)l k ] so if 
7 = DI is a solution to r = 0, then one must have [dLi nt /dj] Dl = P(D)1 = for 
some polynomial P(D).) It must be pointed out that Bias, Deffayet, and Garriga already 
claimed for the existence of solutions of the kind considered in equation ( 125]) in bigravity 
for an arbitrary interaction term between the metrics, where D would be determined 
by the particular interaction term |HD]- The novelty in our study resides in the fact 
that we have concluded that those solutions are solutions of massive gravity which are 
continuous in the parameter space. Moreover, for the particular interaction Lagrangian 
considered in this paper, solution ( 125]) is not the only solution implying an Einstein 
manifold g, as it was the case for the particular theory considered in reference [60J (see 
also reference [55] for a particular model), because any solution of bimetric gravity with 
t^v oc 5^ v would lead to a foreground metric of that kind. 
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3.2.2. Second general solution. The second general solution corresponds to the metric 

dx» dx v = D 2 S 2 {r, t) dt 2 - D 2 N 2 {r, t) dr 2 - D 2 R 2 {r, t) dQ 2 {2) , (29) 
where D and D are two separate constants satisfying 

D - -D + <*(D-l) + 2 

c 3 p-l) + l ' ^ 

and 

-1 - c 3 + c 2 (D - l) 2 - c 3 D 

04 ~ (iJ^iy • (31) 

In this case, the extra term in the equations of motion ( fill behaves as a positive 
contribution to the cosmological constant if c 3 (D — 1) > — 1, leading to a foreground 
universe with accelerating expansion. 

This solution is similar to the previous one, with the only difference that the t-r 
sector of both metrics is conformally related through one constant whereas the angular 
sectors, in which we imposed the symmetry, are related through another constant. This 
fact would not lead to great differences, at least in principle. 

3.2.3. Two special case solutions. There are two additional special-case solutions that 
hold only for a specific relation between the parameters C3 and C4. In particular, we need 
to impose that C4 = —3/ Ac 2 . Both solutions lead to an effect in the modified Einstein 
equations (ITTj) equivalent to a negative contribution to the cosmological constant. For 
the first special case solution, the background metric can be written as 

dx» dx v = D 2 S 2 (r, t) dt 2 - B 2 (r, t) dr 2 - D 2 R 2 (r,t) dtf 2) , (32) 

whereas for the second one 

/ M „ dx" dx v = A 2 (r, t) dt 2 — D 2 N 2 (r, t) dr 2 — D 2 R 2 (r, t) dQ 2 {2) . (33) 

Although these solutions are compatible only with a particular sub-class of models, they 
now allow us to have two possibly different cosmological metrics in the background and 
foreground. 

Regarding the first kind of solutions, where the foreground and background metric 
are given by equations (122"]) and (1521 . respectively, let assume that there are solutions 
with a foreground metric of the FLRW-kind. That is 

g„ v dx^ dx v = a 2 (t) dt 2 - ^ 2 dr 2 - a 2 (t)r 2 dtt 2 2) , (34) 

with t being the so-called FLRW conformal time coordinate (it is not the more usually 
occurring FLRW proper time coordinate) and the scale factor a(t) fulfilling the equation 
of motion ([Til , which, for this symmetry and this coordinate system, can be expressed 
as 

* 2 7 2 

3 Q + A = A — m 2 (D — l) 2 + 8nG Pm . (35) 
a 4 

Thus, although the effect of the graviton mass is equivalent to a negative contribution 
to the cosmological constant, the solution is accelerating if A cff > 0, with A eff = 
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A— m 2 (D— l) 2 . The foreground spacetime (I34p is a solution of massive gravity continuous 
in the parameter space, if we consider a theory with a background metric which can be 



written as (see Appendix C) 



dx" dx v = D 2 a 2 (t) dt 2 - B 2 (r, t) dr 2 — D 2 a 2 (t) r 2 dQ 2 (2) . (36) 

As the function B(r,t) is not constrained by our analysis, it can take any form, not 
necessarily related with a(t). Thus, all the massive gravity theories defined with a 
background metric of the form given by equation fl36l) have FLRW solutions which 
are continuous in the parameter space. One interesting example would be to consider 
B 2 {r, t) = D 2 a 2 {t) / {I — kr 2 ), in such a way that the background metric now has different 
spatial curvature as the physical metric. One could also consider more exotic cosmologies 
compatible with metric fl36l) . 

We can also reverse the logic of the problem and take a massive gravity theory 
defined with an isotropic and homogeneous background metric. That is 

/ M „ dx» dx v = a 2 (t) dt 2 - f^} 2 dr 2 -a 2 (t) r 2 dtt 2 {2) , (37) 

which includes the particular case of massive gravity theories with a de Sitter background 
metric. In this case, requiring the fulfillment of the constraints ( TT8|) (and (I20I) ). it can 
be seen that the foreground metric can be written as 

g, v dx» dx v = ^dt 2 -N 2 (r,t) dr 2 - dQ 2 {2) . (38) 

Note that the background scale factor a(t) does not now fulfill any Friedmann-like 
equation, since the background metric is not constrained but externally specified by 
the theory itself. We can obtain various different kinds of physical cosmologies, being 
described by the foreground metric, by changing the function N(r,t). As already 
mentioned, one particular solution belonging to this class would be that assuming that 
both spaces are FLRW, obtained by fixing N(r, t) = a 2 (t)/[D 2 (l — kr 2 )]. Nevertheless, it 
is probably more interesting to consider the solutions implying an anisotropic expanding 
universe for a massive gravity theory defined by using a FLRW background metric, which 
could even have the maximal symmetry being a de Sitter space. That would be the case 
of any solution obtained by choosing N(r,t) = b 2 (t)/[D 2 (l — kr 2 )], with b(t) ^ a(t). 

On the other hand, those solutions in which the background metric takes the form 
given by metric ( |33l) have similar characteristics as those with ( |32l) . The only difference 
is that in this case one can arbitrarily choose the function appearing in the temporal 
component of the metric. That is, for a theory with a FLRW background metric ([3] 
one has 



a 2 (t) ^ 2 a 2 (t)r 2 



dx» dx v = S 2 (r, t) dt 2 - —^-L— dr 2 - — ^ dtt 2 {2)} (39) 

which again can be fixed to be homogeneous and isotropic but it also allows the 
consideration of more exotic situations. 

Finally, it must be pointed out that the solutions presented here are, on one hand, 
more general that those already studied in the literature (for instance they can describe 
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anisotropic cosmologies in one sector) and, on the other hand, more restrictive than 
other solutions, as they are only solutions for a particular model (with the parameters 
appearing in the Lagrangian fulfilling a particular relation). 

3. 3. Limit procedure without vanishing background matter 

Let us briefly consider a slightly different limit which retains some effects of the 
background matter, that is k — > with e nonzero and fixed. In this case the action 
becomes 

S = - J d A x^g~ {R(g) + 2 A - 2m 2 L int } + S {m) + eS (m) , (40) 

with dependence on the background metric f^ u both in the obvious place L int (f,g) and 
in the action e S( m ) for the background matter fields, which now have only an extremely 
indirect influence on the foreground physical sector. In this case the equation of motion 
for the background metric f^ u becomes a purely algebraic one 

T^ = - e -^f(™K, (41) 
mr 

in terms of the background matter. One still has some equation constraining the 
background metric, and there will be some constraint intertwining the background and 
foreground metrics, and this is basically unavoidable in any reasonable limit of bimetric 
gravity that is based on tuning the parameters in the bigravity action to specific values. 
Nevertheless, it should be noted that one would not be able to directly detect the 
existence of any matter content in the background. Therefore, although the theory 
resulting when considering the limit k — > (with e nonzero and fixed) has a very different 
motivation that massive gravity, it seems that we could not distinguish between them by 
measuring their physical "foreground" consequences. In other words, we cannot think 
in any physical prediction which would be affected by considering that f^ u is given by 
the theory, or constrained by some matter invisible to us. 



4. Summary and Discussion 

In this work we have explored the cosmological solutions for a massive gravity theory 
when viewed as a limit of taking a vanishing kinetic term in a bimetric theory, paying 
particular attention in the way this limit is taken. We have used a vierbein formulation 
based on that developed by Volkov in reference [40] (see also references [621 [36J EJ [19] ) , 
a formalism that proved to be very powerful in the treatment of this type of calculation. 

A first step is to realize that the solutions of massive gravity, taken as the limit of a 
bigravity theory, are in general more constrained than the general solutions of massive 
gravity. That is because in massive gravity there is only one set of equations, whereas 
if one considers this theory as a particular limit of bimetric gravity, then additional 
constraints must be taken into account. Thus, one cannot recover complete continuity 
in the physical predictions of the theory, since the solutions continuous in the parameter 
space of the theory are not the complete class of solutions. 
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We have argued that massive gravity can be recovered from bimetric gravity by 
suitably taking the limit for k — > in the action describing the latter. We have shown 
that this limit implies that the modification to the equation of motion introduced by a 
non-zero mass of the graviton is equivalent to introducing an extra contribution to the 
cosmological constant which, at this stage, can be positive or negative. An important 
fact to stress is that no restriction on curvature or dynamics of the background metric 
has been imposed. 

In particular, we have focused on ghost-free cosmologies with both metrics, 
foreground and background, being spherically symmetric. In the first place, we have 
shown that there are no non commonly diagonal solutions continuous in the parameter 
space. In the second place, we have found two kinds of solutions for any ghost-free theory. 
The first kind of solutions implies that the two metrics are proportional to each other, 
with the graviton mass interaction Lagrangian producing an extra contribution to the 
cosmological constant that can be either positive or negative, depending on some relation 
between the parameters of the theory. The net result is that this solution describes an 
accelerating or decelerating universe which has the same symmetry as the background 
metric. The second kind of solutions implies a more complicated relation between the 
two metrics, but we still obtain a positive cosmological constant for certain values of 
C3. Moreover, in the third place, we have also seen that two more kind of solutions can 
be obtained by considering a very specific relation between the parameters C3 and C4, 
and that these solutions are equivalent to inducing an extra negative contribution to 
the cosmological constant in the modified Einstein theory. These solutions allow us to 
consider different cosmological metrics, for the background and foreground, related by 
some unconstrained arbitrary functions. In the first kind of solutions, when considering 
a spherically symmetric background metric and a FLRW physical metric, with a scale 
factor fulfilling the Einstein equations, we can obtain a different spatial curvature for the 
two metrics, by particularly tuning the function appearing in the radial component of 
the metric. If instead we consider the background metric to be FLRW, the (background) 
scale factor is now not forced to fulfill any Friedmann equations since this metric does 
not have a dynamics. A particular interesting case is to consider a massive gravity theory 
with an isotropic and homogeneous background metric which can lead to anisotropic 
expansion of the foreground space. On the other hand, in the second kind of solutions, 
we can obtain different cosmologies by changing the function in the temporal component 
of the metric. 

It must be emphasized that the solutions obtained in this paper are solutions of 
massive gravity continuous in the parameter space of this theory; that is, they are 
solutions simultaneously of both massive gravity and bimetric gravity in the limit of 
a vanishing kinetic term. Thus, these solutions are completely different from those of 
massive cosmologies presented in references j35j [36j E3 EHJ HI] , which are only solutions 
of massive gravity, and from those bimetric cosmologies studied in references [4T H I4TI |4~2] , 
which are not solutions when the limit of vanishing kinetic term is taken. On the 
other hand, these solutions must also not be confused with the solutions obtained in 
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references [4"2"l |4"0"] where a different limit of bimetric gravity is taken to recover massive 
gravity, namely k — > oo. 

On the other hand, in view of the mentioned results, one could consider that the 
consideration of massive gravity as a limit of bimetric gravity leads to a theory which is 
equivalent to general relativity. Nevertheless, if this conclusion remain unchanged by the 
consideration of perturbative effects in the metric (note that the effective cosmological 
constant appears when considering the particular solutions and not directly in the 
action), then that would not be a problem of our treatment but of the theory of massive 
gravity itself. In that case one should conclude that: either (i) massive gravity is 
continuous on the parameter space only when it is equivalent to general relativity (with a 
particular value for the cosmological constant), or (ii) e cannot be zero and the dynamics 
of our universe is affected by some background invisible matter, or (iii) massive gravity 
cannot be seen as a particular limit of bimetric gravity. We have briefly considered 
option (ii) in this paper, but there is a rich phenomenology of models that might be 
explored in this case. It could be interesting to consider whether the kinds of models that 
should be considered could be constrained by some property of the hidden matter sector. 
For example, could it be reasonable to require that the classical energy conditions be 
satisfied in this sector? Should there be some relationship between the standard model 
of particle physics in both sectors, or can they be completely independent? Or, could 
this hidden matter be the "mirror matter" speculated to restore complete symmetry in 
the fundamental interactions [63] [M]? 

In any case, in many ways it seems that bimetric theories and massive gravity 
should be kept in conceptually separate compartments. While there is no doubt that 
the relevant actions are very closely related, as we have seen treating massive gravity 
as a limit of bimetric gravity is fraught with considerable difficulty In fact, whereas in 
massive gravity L int simply gives mass to the graviton, in bimetric (and multimetric) 
gravity this term is actually describing an interaction between two (or more) metrics. 
Therefore, bimetric gravity could be considered as a model of a bi-universe, in which two 
different physical worlds are coexisting, and interacting only through the gravitational 
sector (at a kinematical level the "analogue spacetime" programme can be used to 
provide examples of multi- metric, though not multi-gravity, universes [65 | l66 l l68 | 167]). 
Moreover, such an interpretation would lead to a new multiversal framework when 
considering multimetric gravitational theories [19], which would be very different to 
those previously considered in the literature, since it is not resulting from quantum 
gravity [69J, inflationary theory [70], string theory [71], or general relativity [72] , but it 
is a consequence of abandoning Einstein's theory. This type of multiverse is also rather 
different from the more exotic multiverse concept developed in references [73] 174"] . 
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Appendix A. Some identities 

In this appendix we collect some purely algebraic results. For a general n x n matrix X 
the symmetric polynomials ei(X) are defined by 

n, 

J^A* e;(X) = det(I + AX). (A.l) 

i=0 

Thus, the symmetric polynomials can be recovered iteratively from Newton's identity 

= ) Ec- 1 )' -1 e ^ x ) fc^i. ( A - 2 ) 

3=1 

taking into account eo(x) = 1. Since 

, /it ^ u detpT + AI) A n det(I + X^X) . . , 

de t (I + AX-) = ^ U (A.3) 



we have 



—in e n _i(X) 



which is a purely algebraic result which allows us to rewrite the interaction term 
in bimetric gravity in various useful ways [17] (see Appendix B for more details). 
Furthermore, since 

n n 

^e^I + = det((l + A)I + eAX) = + A) n " i (eA) i e;(X) (A.5) 

i=0 i=0 

n n—i / \ n n / .\ 

= EEiVr^^EEu-O^w (A - 6) 

8=0 j=0 ^ J ' 8=0 fc=i ^ 7 

n k / . \ n * / 7 \ 

k=0 8=0 V 7 8=0 fc=0 v 7 

we have the "shifting theorem" 

ei p + do = e (" : •) <* * w = e (": J) «* ^< a ')- < A - 8 ) 

k=0 ^ 7 fc=0 ^ 7 



Appendix B. The interaction term in bimetric gravity 

In this paper we have used an expression for the graviton mass term, equation (jSJ), which 
emphasized the fact that ghost-free massive gravity only includes three parameters more 
than general relativity. The consideration of such a L mt in massive gravity leads to the 
same equations of motion as those that can be obtained using a formulation of the 
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interaction term in terms of 7 as used in reference [5]. In this appendix we will show 
that this is also the case in bimetric gravity when one cannot throw away the term 
depending only on the background metric. That is, both formulations are equivalent 
also in bimetric gravity, at least when considering both a foreground and background 
cosmo logical constant. 

We note that the ghost-free foreground-background interaction terms can in all 
generality be written as 

4 4 

V^g L int = ki e ^ 7 ) = v 7 -/ e ^~ 1 )- ( B - 1 ) 

i=0 i=0 

Here we have used the explicit algebraic symmetry between ei(X) and e n -i(X~ l ) 
to exhibit an explicit interchange symmetry between foreground and background 
geometries. (See reference [32].) Furthermore, in view of the fact that 7 = I — K, 
the shifting theorem yields 

e 1 ( 7 )=4-e 1 (K); 

e 2 ( 1 ) = 6-3e 1 (K)+e 2 (K); 

e 3 ( 7 ) = 4 - 3 ei (K) + 2e 2 (K) - e 3 (K); 

e 4 ( 7 ) = l- ei (K) + e 2 (K) - e 3 (K) + e 4 (K). 

Consider the original "minimal" Lagrangian for generating a graviton mass [13] : 



J old minimal 



tr 



(v 7 ^ 1 /) " 3 = tr( 7 ) - 3 = ei ( 7 ) -3 = 1- ei (K). (B.2) 

Using the last of the shifting theorem equivalences from Appendix A , which relates 64(7) 
with the polynomials in K, we see 

^oidminimai = 64(7) - ^ 2 (K) + e 3 (K) - e 4 (K) . (B.3) 

When written in this way the 64(7) term corresponds in bimetric gravity to a background 
cosmological constant, and in massive gravity to an irrelevant constant. The e 2 (K), 
e 3 (K), and e 4 (K) terms are manifestly quadratic, cubic, and quartic in K. So insofar 
as one is only interested in giving the graviton a mass the quantity 

Aicwminimal = _ ^(K), (B-4) 

does just as good a job (in fact arguably a better job) than the original minimal mass 
term. This argument generalizes, we can use the shifting theorem to rewrite the general 
interaction term as 

4 4 

Lint = J2 ki = J2 ki e * ( R5 ) 

i=0 i=0 

Explicitly separating off the first two terms and using again the last of the shifting 
theorem equivalences we see 

4 

Lint = h + fci {1 - e 4 (7) + e 2 (K) - e 3 (K) + e 4 (K)} + ^ k e^K). (B.6) 

i=2 
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So we have 

L^t = (k + k) e {K) + (k 2 + k) e 2 (K) + (k 3 - k) e 3 {K) + {~h + k) e^K) 

-ke^(i). (B.7) 

This eliminates (or rather redistributes) the e\(K) term. For current purposes it is 
now useful to split off the top and bottom terms (corresponding to foreground and 
background cosmological constants) and deal with them separately. We have explicitly 
absorbed them into the bigravity kinetic terms. After factoring out an explicit m 2 , this 
finally leaves us with the three-term interaction Lagrangian we have used in the paper, 
equation (|5]). This is 

L int = e 2 {K) - c 3 e 3 {K) - c 4 e A {K). (B.8) 

Thus, we have shown explicitly that the interaction term (IB.8j) describes the same 



bimetric gravity theory as that considered in ( IB. II) when one takes into account a 
cosmological constant for each metric. 

On the other hand, there is still a hidden "symmetry" between foreground and 
background. We can define an equivalent Kf for the background metric with respect to 
the foreground metric via 

Kf — I — 7 _1 = I + (I — K)- 1 = -K(l -K)~\ (B.9) 

and equivalent ly 

K = -K f {l-K f y\ (B.10) 

Note that 

e 2 (K) =e 2 (K f ) + 0(K 3 f ). (B.ll) 

Therefore, as could have been suspected from equation fIB.ll) . one can equivalently 
consider that the interaction term is giving mass either to the graviton related with f^ v 
or to that of g^. In fact, the ghost-free bimetric gravity is giving us 7 degrees of freedom 
(14 considering also the conjugate momenta) to distribute between both gravitons (see 
reference [5]), without any particular preference as to which. Here resides the great 
conceptual difference between massive gravity and bimetric gravity, since whereas in 
the first theory L mt is a term whose purpose is giving mass to the graviton (the only 
graviton present in this theory), in bimetric gravity the non-existence of a preferred 
metric leads one to interpret L; nt merely as an interaction term. In fact, one could even 
think of some kind of democratic principle for bimetric theories, by interpreting the 
degrees of freedom to be distributed in such a way that we have two massless gravitons, 
and an interaction between the two metrics mediated by one vectorial and one scalar 
degree of freedom. 
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Appendix C. Spherically symmetric solutions 



Taking into account the spherically symmetric metrics (1221) and ( 1231) in equation ( TT71) . 
the non-vanishing components of r^ v can be written as [ID] 



A B C 2 
Jn + N 2 



3A 



+C 4 1 - 





2AU 
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In particular this implies 



T 



BS-AN 
~CN 



{C.l) 



(C.2) 



which greatly simplifies some calculations. We must find solutions for all 
components being set equal to zero. Let us start by requiring that r° r = 0. 
can be obtained in either of two ways: 

(i) C^O, with U(r,t) = D R(r,t), and an appropriate constraint on D. 

(ii) C = 0. 



(C.3) 



(C.4) 



(C.5) 

these 
This 



Appendix C.l. Non-diagonal background metric 

If we wish to consider a non-diagonal background metric, then we must require 
U(r,t) = D ■ R(r, t) in order to have r° r = 0, where D is a constant such that 
3 - 2D - c 3 (3 - AD + D 2 ) 



c 4 



(C.6) 



(D-l) 2 

and, of course, D ^ 1. Substituting this into r r r — r% we find that r r r — r 6 ' q vanishes 
when 

D — 2 



D-l 



(C.7) 
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c 4 {D _ 1)2 ■ (C.8) 

The components of t m „ then reduce to 

tV = D(D + ^] (5V (C.9) 



^N 2 NS / 

Thus, one must require A = —C 2 S/(NB) to have a non-trivial situation. The 
interaction term now reads 

Lm = D(D- 1) (d| + £ - I) = - 1). (C.IO) 

Replacing the relation A = —C 2 S/(NB) in equation f|2^|) for the background metric, 
one has 



dx" dx" = — - ij dt 2 + 2— (-- 

C 2 

'|2 T1 J U ' _J " " U "(2) 5 



52 -^5 + 1 dr 2 - D 2 r 2 dn 2 (2) , (en; 



dx^ dx* = — - 1 — dt - 5dr — -D 2 i? 2 dfiL, (C.12) 



which is non-Lorentzian. In fact this metric can be written as 

C 2 \ /CS* \ 2 
_-lj^_df-Sdr) — «^ (a) , 

which manifestly has unphysical null signature (0, — sign [i? 2 — C 2 ], —1,-1). Therefore 
there is no physical solution of this kind. 

Appendix C.2. Diagonal background metric 

We now set C = 0. In order to have r r r — r° = we have two possibilities, either 
U(r, t) = D ■ R(r, t) with D ^ 1 and such that 

3-2D-c 3 (3-AD + D 2 ) 
c 4 = (D-l) 2 ' ( } 

as before ("case I"), or B S = A N ("case II"). 

Case I: Solutions for particular models. If we consider the first case, U = D R, then 
in order to have also t£ — Tq = 0, there are three options: 

(i) C3 takes the same value as in the previous subsection (when considering a 
non-diagonal background metric): this implies the same consequences. The 
background metric is non-Lorentzian now specifically with unphysical null signature 
(0,-1,-1,-1). 

(ii) A = D S: In this case we have 

t ^ = B(2 + c 3 (D-1))(D-1)D s ^ (C14) 
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thus, the theory should have 

2 3 
c 3 = [, and c 4 = - , (C.15) 

to have solutions. The Lagrangian is L int — (D — l) 2 , which fulfills the Bianchi- 
inspired constraint, 
(iii) B = DN: Now 

^ = Ai2 + C3[ D- m D-l)D ^ (c 16) 

which must vanish. Therefore, we have again c 3 and C4 given by equation ( 1C15|) . 
and the same L int as in the previous case. 

Case II: Solutions for all the ghost-free models. We now consider A = B S/N. It can 
be seen that there are two cases in which r r r — r e g = 0. These are: 

(i) B = N U/R: We then have 

(g - u)u[(3 - 3c 3 - c,)R 2 + (3c 3 + 2c,)R U - c,U 2 } 
T v = ^4 <>V (C.17) 

These quantities vanish for U = D ■ R, with D such that 

C = 5±§^, (C,8) 

which implies L int = (3 + c 3 (-D — 1))(D — l) 2 . Thus, we should require D ^ 1 to 
have a (non-trivial) massive gravity theory. 

(ii) Consider 

3-3c 3 -c 4 + (2c 3 + c 4 - 1 ) t //J i 

1 - 2c 3 - c 4 + (c 3 + 04)^//? V ; 

In this case we have 

„ U_ [(3-3c 3 -c A )R-(l-2c 3 -c 4 )U] 
Tu R* [(2c 3 + C4-l)-f/(c 3 + c 4 ) 2 

x \(cl - c 3 + c 4 + l)i? 2 + (c 3 - 2c 2 . - 2c±)RU + (c 2 , + c 4 )f/ 2 ] 5^, (C.20) 

which again vanishes if U = D ■ R, but now we must have .D such that 

c§(D-l) 2 + C3 (D-l) + l 
c 4 - • (C.21) 



The mass term now reads 

int " c 3 p-i) + i- 



(C.22) 



It should be noted that there is yet one more formal solution for the equations 
r^ u = 0, but one which at the end of the day implies B = 0. Thus, for that solution 
the background metric is again non-Lorentzian and unphysical, this time with null 
signature (+1, 0, —1, —1). 

This completes our consideration of the various explicit constraints on the background 
geometry arising from the equation r = in the case of spherical symmetry. 
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